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Abstract 

Published  accounts  of  abuse  testing  of  lithium-ion  cells  and  components  are  summarized,  including  modeling  work.  From  this  summary,  a 
set  of  exothermic  reactions  is  selected  with  corresponding  estimates  of  heats  of  reaction.  Using  this  set  of  reactions,  along  with  estimated 
kinetic  parameters  and  designs  for  high-rate  batteries,  models  for  the  abuse  behavior  (oven,  short-circuit,  overcharge,  nail,  crush)  are 
developed.  Finally,  the  models  are  used  to  determine  that  fluorinated  binder  plays  a  relatively  unimportant  role  in  thermal  runaway. 
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1.  Introduction 

With  increasing  interest  in  lithium-ion  batteries  for  auto¬ 
motive  applications,  there  is  a  need  to  better  understand  the 
abuse  tolerance  of  these  batteries.  For  example,  that  fluori¬ 
nated  binders,  such  as  polyvinylidene  fluoride  (PVDF),  react 
exothermically  with  lithium  is  of  some  interest  because  use 
of  fluorinated  binders  is  desirable.  Fluorinated  binders  per¬ 
form  well  and  are  easy  to  use.  So,  the  question  naturally 
arises,  do  fluorinated  binders  contribute  to  thermal  runaway? 
Or,  because  fluorinated  binders  react  in  the  negative  elec¬ 
trode  only  above  200  °C,  is  the  reaction  a  consequence  of 
thermal  runaway?  If  so,  to  what  extent  does  it  contribute  to 
the  overall  heat  release  and  thus  aggravate  thermal  runway? 

The  question  as  to  the  role  of  fluorinated  binders  in 
thermal  runaway  can  be  addressed  through  mathematical 
modeling.  Through  modeling,  the  various  reactions  contri¬ 
buting  to  thermal  runaway  can  be  decoupled  and  so  clarify 
what  reactions  cause  thermal  runaway.  To  achieve  this  end, 
the  model  should  include  the  chemical  reactions  and  the 
mechanisms  for  heat  conduction  and  dissipation.  Modeling 
is  practically  indispensable  in  resolving  this  issue,  since  it 
enables  a  broad  range  of  compositions  and  conditions  to  be 
explored,  without  the  need  for  assembling  and  testing  an 
extensive  number  of  large  batteries. 

Here,  simulation  results,  based  on  a  mathematical  model 
for  abuse  tolerance  (oven,  short-circuit,  overcharge,  nail  and 
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crush),  are  presented.  The  simulations  are  based  on  the 
thermal  behavior  of  lithium-ion  battery  materials.  Before 
presenting  the  model,  reviews  of  the  abuse  behavior  of 
lithium-ion  batteries  and  modeling  work  are  presented. 

2.  Survey  of  abuse  behavior  of  lithium-ion  batteries 

The  following  abuse  tests  are  widely  used  to  characterize 
the  abuse  tolerance  of  lithium-ion  cells. 

(a)  Oven  test:  This  test  simply  involves  exposing  the 
battery,  at  some  initial  temperature,  to  a  higher 
temperature.  For  consumer  batteries,  an  oven  tempera¬ 
ture  of  150  °C  is  used  (UL2054). 

(b)  Short-circuit:  A  low  resistance  (<5  mQ)  is  connected 
across  the  terminals  of  the  battery.  The  battery  may  be 
preheated.  In  this  test,  current  flows  through  the  battery 
generating  heat.  The  battery  is  heated  internally  due  to 
current  flow,  but  the  external  circuit  can  dissipate  heat 
also. 

(c)  Overcharge  test:  Current  is  forced  through  the  cell  up  to 
some  limiting  voltage.  Heat  is  generated  by  electro¬ 
chemical  reactions  and  by  current  flowing  through  the 
cell. 

(d)  Nail:  A  nail  is  forced  through  the  battery  at  a  prescribed 
rate  (such  as  8  cm/s).  Heat  is  generated  by  current 
flowing  through  the  cell,  and  by  current  flowing 
through  the  nail.  Initially  the  nail  is  positioned  outside 
of  the  battery  wall  and,  when  the  test  begins,  is  forced 
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Nomenclature 

<2s?n  specific  area  of  negative  active  (cm 2/g) 

Bmw  molecular  weight  of  binder  (g/mol) 

Li  site  concentration  in  negative  (mol/cm3) 

Cq  moles  of  oxidizable  solvent  per  unit  volume  of 
electrolyte  (mol/cm3) 

cR  moles  of  reducible  solvent  per  unit  volume  of 
electrolyte  (mol/cm3) 

cs  moles  of  solvent  per  unit  volume  of  electrolyte 
(mol/cm3) 

Eneg  voltage  of  negative  versus  lithium  (V) 

£pOS  voltage  of  positive  versus  lithium  (V) 

F  Faraday’s  constant  (C/eq.) 

Hub  enthalpy  of  lithium/binder  reaction  (J/mol) 

Hus  enthalpy  of  lithium/solvent  reaction  (J/mol) 
HNB  enthalpy  of  negative/binder  reaction  (J/mol) 
Hns  enthalpy  of  negative/solvent  reaction  (J/mol) 

HP  enthalpy  of  positive  decomposition  reaction 
(J/mol) 

Hs  heat  of  reaction  for  solvent  decomposition  (J/mol) 
Hsu  (J/mol) 

i  current  density  per  area  of  electrode  (A/cm2) 

iu  current  density  corresponding  to  lithium  metal 

deposition  (A/cm2) 

ioc  overcharge  current  density  (A/cm2) 

I  inert  material  (mol) 

kd  rate  constant  for  lithium  deintercalation  under 

short  circuit  conditions  (mol/cm2) 
ku  scaling  factor  to  relate  moles  of  Li  to  area 
(cm2/mol) 

&EiB  rate  constant  for  Li/binder  reaction  (s-1) 
kus  rate  constant  for  Li/solvent  reaction  (s-1) 
kmN  mass-transfer  coefficient  for  lithium  diffusion  in 
negative  (cm/s) 

kNB  rate  constant  for  negative/binder  reaction  (s-1) 

kNS  rate  constant  for  SEI  formation  (cm3/(mol  s)) 

kp  positive  decomposition  rate  constant  (s-1) 

ks  solvent  decomposition  rate  constant  (s-1) 

&sei  rate  constant  for  SEI  decomposition  (s-1) 

/N  coating  weight  per  area  of  negative  (g/cm2) 

Le  thickness  of  SEI  layer  (cm) 

LE  minimum  thickness  of  SEI  layer  (cm) 

Lt  thickness,  i  =  N,  P,  or  sep  (cm) 

Ln  thickness  of  the  negative  coating  (cm) 

Lp  thickness  of  positive  coating  (cm) 

Lpos  thickness  of  positive  (=LP)  (cm) 

Lsep  thickness  of  separator 

mE  moles  of  SEI  per  unit  area  of  negative 
(mol/cm2) 

mu  moles  of  Li  metal  per  unit  area  of  negative 
(mol/cm2) 

mN  mass  of  the  negative  active  material  (g/cm2) 
raPO  moles  of  positive  active  material  in  the  oxidized 
state  per  unit  area  of  positive  (mol/cm2) 


mPR  moles  of  positive  active  material  in  the 
reduced  state  (mol/cm2) 

^crushed  number  of  cells  “crushed”  in  crush  test 
Ashorted  number  of  cells  shorted  by  nail  in  nail  test 
An  mol  of  Li  per  cm2  of  cell  area  (mol/(cm2  s)) 

Rub  rate  of  lithium/binder  reaction  (mol/(cm2  s)) 

Rus  rate  of  lithium/solvent  reaction  (mol/(cm2  s)) 

Rnb  rate  of  negative/binder  reaction  (mol/(cm3  s)) 

Rns  rate  of  negative/solvent  reaction  (mol/(cm2  s)) 
Rp  rate  of  positive  decomposition  (mol/(cm2  s)) 

Rs  rate  of  solvent  decomposition  (mol/(cm3  s)) 
Rsei  rate  of  SEI  decomposition  (mol/(cm2  s)) 

Rsep  separator  resistance  (Q  cm2) 

SN  heat  source  term  for  negative  (W/cm2) 

Sp  heat  source  term  for  positive  (W/cm2) 

Ss  heat  source  term  for  separator  (W/cm2) 

wB  weight  fraction  of  binder  in  negative 
wN  weight  fraction  of  negative  active  material 
in  negative  coating 

xh  fraction  of  carbon  sites  filled  in  bulk 

vs  fraction  of  carbon  sites  filled  with  Li  at  surface 

Greek  letters 

eN  porosity  of  negative 

Sp  porosity  of  positive 

vN  volume  fraction  of  negative  active  in 
negative  coating 

Dps  stoichiometry  for  solvent/oxygen  reaction 

us  stoichiometry  for  electrochemical  oxidation 

of  solvent 

Dus  stoichiometry  for  lithium/solvent  reaction 

Dub  stoichiometry  for  lithium/binder  reaction 

uNB  stoichiometry  for  negative/binder  reaction 

uNS  stoichiometry  for  negative/solvent  reaction 

pB  molar  density  of  SEI  (mol/cm3) 
ot  conductivity  of  layer  or  material 
(Q-1  cm-1) 


through  the  battery  wall  and  into  the  battery  at  a 
constant  speed.  As  the  nail  moves  forward,  forming 
direct  shorts  between  adjacent  electrode  pairs,  the 
current  flowing  through  the  nail  itself  decreases. 

(e)  Crush  test:  A  bar  is  used  to  press  down  on  the  battery 
until  a  short-circuit  initiates.  Heat  is  generated 
throughout  the  battery  due  to  discharge,  and  locally 
through  the  electrode  pairs  that  are  physically  shorted. 

The  overcharge  test  can  be  the  most  severe  as  additional 
energy  is  added  to  the  cell.  The  nail  and  crush  tests  can  be 
enhanced  by  using  overcharged  batteries.  The  nail,  crush, 
and  short-circuit  tests  are  somewhat  similar,  but  have  sig¬ 
nificant  differences.  In  these  three  tests,  heat  is  generated  by 
rapid  discharge  of  the  battery.  The  short-circuit  test  allows 
the  cell  to  discharge  uniformly  throughout  the  volume  of  the 
cell.  The  nail  test  allows  the  cell  to  discharge  uniformly,  but 
focuses  heat  (due  to  the  nail)  locally.  In  addition,  the  nail 
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test,  since  the  nail  penetrates  the  casing  of  the  battery,  allows 
volatile  electrolyte  to  leave  and  react  with  air.  The  crush  test 
allows  the  cell  to  discharge  uniformly  from  the  short,  but 
heat  is  generated  locally  in  the  cell  due  to  the  short.  If  the 
battery  is  designed  to  form  a  good  internal  electrical  contact 
when  crushed,  then  the  crush  test  will  behave  like  a  short- 
circuit  test. 

Only  a  few  reports  on  the  abuse  behavior  of  large  cells  are 
available.  Most  of  the  published  accounts  of  abuse  behavior 
refer  to  small  (<1.5  Ah)  lithium-ion  cells. 

Fouchard  et  al.  [14]  reported  the  response  of  AA-size, 
lithium-ion  cells  to  heating  in  an  oven.  They  found  the  onset 
temperature  for  thermal  runaway  varied  inversely  with  the 
degree  of  lithiation  of  the  negative  electrode  (the  greater  the 
degree  of  lithiation  the  lower  the  onset  temperature),  while 
the  thermal  stability  of  the  cell  increased  with  increasing 
lithiation  of  the  cathode.  They  reported  that  the  composition 
of  the  electrolyte  affected  the  thermal  stability  of  the  cell 
(ethylene  carbonate/dimethyl  carbonate  (EC/DMC)  was 
more  stable  than  PC/DME),  and  that  LiMn204  cathodes 
offered  better  thermal  stability  than  LiCo02  or  LiNi02 
cathodes. 

Kitoh  and  Nemoto  [23]  report  the  results  of  external  short- 
circuit,  overcharge,  nail  penetration,  and  oven  tests  with 
cylindrical,  100  Wh,  LiMn204/carbon  batteries.  The  bat¬ 
teries  contained  shutdown  separators,  fuses,  and  safety 
vents.  In  a  130  °C  oven  test,  the  battery  can  temperature 
seemed  to  stabilize  near  the  oven  temperature,  but  after 
about  80  min,  the  can  temperature  jumped  to  ~190  °C.  At 
this  point,  the  battery  vented,  but  did  not  ignite.  About 
40  min  into  the  test,  the  battery  exhibited  a  large  increase 
in  internal  resistance,  most  likely  due  to  melting  of  the 
separator.  In  a  short-circuit  test,  an  external  load  of  6.0  mQ 
was  used  compared  to  a  battery  impedance  of  ^5  mQ.  Once 
shorted,  the  battery  discharged  at  ~15C  rate  for  ~80  s;  this 
corresponds  to  ~33%  of  the  total  capacity.  At  ~80  s,  the  cell 
vented  and  the  battery  continued  to  discharge  at  a  decreasing 
rate.  Most  likely,  the  decrease  in  current  was  partly  due  to 
shutdown  of  the  separator.  If  the  short-circuit  test  was 
carried  out  with  an  external  load  of  0.2  mQ,  the  current 
was  immediately  cut-off  due  to  a  fuse  inside  the  battery.  In 
an  overcharge  test,  a  1 C  current  was  applied  with  a  voltage 
limit  of  10  V.  When  the  can  temperature  reached  ~100  °C, 
the  current  dropped  dramatically  due  to  shutdown  of  separa¬ 
tor,  and  the  test  was  discontinued.  In  a  nail  penetration  test 
(nail  speed  1  mm/s),  the  cell  vented  immediately  after  the 
nail  was  inserted.  The  battery  reached  a  maximum  tempera¬ 
ture  of  380  °C,  but  did  not  ignite. 

Leising  et  al.  [26]  examined  the  overcharge  and  short- 
circuit  behavior  of  1.5  Ah  prismatic  lithium-ion  cells.  With  a 
thermocouple  placed  inside  the  cell,  they  measured  the 
internal  temperature  versus  time  and  found  it  increased  to 
a  maximum  of  ~132  °C  (melting  point  of  polyethylene 
separator)  whereas  the  outside  surface  of  the  case  reached 
a  maximum  temperature  of  94  °C.  In  overcharge  tests,  they 
confirmed  the  finding  of  Tobishima  and  Yamaki  [24]  that 


low  charging  rates  (<C/5)  did  not  result  in  thermal  runaway, 
while  higher  rates  (C)  did.  By  constructing  cells  with 
varying  ratios  of  cathode  to  anode,  they  found  the  amount 
of  charge  leading  to  thermal  runaway  depended  only  on  the 
amount  of  cathode  material;  runaway  occurs  near  complete 
delithiation  of  the  cobalt  oxide.  They  found  that  cell  rupture 
occurred  when  the  internal  cell  temperature  exceeded 
~190  °C  and  suggest  that  rupture  is  due  to  melting  of 
metallic  lithium  formed  at  the  negative  electrode.  However, 
if  the  cell  is  designed  so  lithium  metal  cannot  deposit,  it  will 
still  rupture,  but  at  a  higher  temperature. 

In  overcharge  studies  of  18,650-size  cells  by  five  different 
manufacturers,  Saito  et  al.  [25]  found  that  the  heat  released 
on  overcharge  was  almost  proportional  to  the  charging 
current.  They  found  that  with  a  0.2C  charging  rate,  at  around 
160-180%  overcharge,  the  heat  generated  in  cells  with 
LiCo02  cathodes  increased  dramatically;  a  cell  with  an 
LiMn204  cathode  exhibited  a  big  increase  in  heat  generation 
at  about  110%  overcharge.  With  further  overcharging,  the 
heat  generated  in  the  cell  became  equal  to  the  electrical 
energy  input.  Using  a  square  wave  to  overcharge  the  cell, 
they  showed  that  shutting  off  the  current  flow  to  the  cell 
stopped  heat  generation.  They  suggest  that  thermal  runaway 
is  the  result  of  solvent  oxidation  products,  formed  above 
4.6  V,  that  react  with  the  positive  electrode  material. 

Maleki  et  al.  [27]  looked  at  the  overcharge  behavior  of 
different  lithium-ion  cells  with  the  objective  of  elucidating 
the  roles  of  the  positive  and  negative  electrodes.  Accelerat¬ 
ing  rate  calorimetry  (ARC)  studies  of  the  full  cells  (type  A) 
and  their  individual  electrodes  indicate  the  positive  electrode 
is  responsible  for  thermal  runaway  at  50  and  100%  state  of 
charge  (SOC),  while  the  negative  electrode  is  responsible  for 
thermal  runaway  at  150  and  200%  SOC.  Above  100%  SOC, 
the  onset  temperature  of  chemical  reactions  for  the  negative 
electrode  decreased  from  >180  to  ~80  °C.  This  suggests  that 
metallic  lithium  plays  an  important  role  in  thermal  runaway 
during  overcharge. 

Biensan  et  al.  [4]  report  nail  tests  done  on  4/5  A-size  cells 
with  different  negative  electrode  binders  and  positive  elec¬ 
trode  materials.  They  found  the  safe  voltage  was  higher 
with  cobalt  oxide  cathodes  than  with  nickel  oxide  cathodes, 
and  with  non-fluorinated  binder  than  with  polyvinylidene 
fluoride  binder. 

Dahn  [34],  using  a  nail  with  an  embedded  thermocouple, 
found  that,  with  18,650-size  LiCo02  cells,  the  nail  tempera¬ 
ture  exceeds  600  °C  when  forced  slowly  into  the  cell  to  a 
depth  <4.5  mm.  However,  with  a  fast,  deep  nail  penetration 
(7.5  mm),  the  nail  temperature  did  not  exceed  140  °C. 

Tobishima  and  Yamaki  [24]  explain  that  the  abuse  toler¬ 
ance  of  lithium-ion  cells  depends  on  the  relative  rates  of  heat 
removal  and  heat  generation.  When  a  cell  cannot  transfer 
heat  to  its  environment  at  a  rate  equal  to  or  higher  than  the 
rate  of  heat  generation,  it  will  undergo  thermal  runaway. 
They  also  present  results  from  various  abuse  tests  of  small 
lithium-ion  cells.  In  the  overcharge  test  of  prismatic  cells  (Al 
case,  600  mAh),  cells  did  not  vent  at  1  or  1.5C  rate  of 
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overcharge,  but  caught  fire  at  a  2 C  rate  overcharge.  In  oven 
tests  of  cylindrical  cells  (1270  mAh),  the  cells  were  stable  at 
150  °C  but  smoked  at  155  °C.  Nail  penetration  tests  with 
prismatic  cells  (835  mAh)  did  not  overheat;  however,  cells 
that  were  slightly  overcharged  (0.03  V  higher  than  standard 
voltage)  did  overheat  and  give  off  smoke.  Similarly,  a  crush 
test  of  a  cylindrical  cell  (720  mAh)  charged  according  to 
the  manufacturer’s  recommendation  did  not  smoke,  but  an 
overcharged  cell  did  catch  fire.  The  fire  is  attributed  by  the 
authors  to  lithium  metal  formation. 

While  there  is  little  information  published  about  abuse 
testing  of  lithium-ion  cells,  there  is  virtually  no  material 
available  on  abuse  testing  of  modules.  The  concern  is  that 
thermal  runaway  of  a  single  cell  will  induce  runaway  of 
an  entire  module.  Modules  envisioned  for  hybrid  electric 
vehicles  (HEVs)  consist  of  at  least  10  cells,  so  runaway  of  a 
module  is  potentially  at  least  an  order  of  magnitude  more 
energetic  than  thermal  runaway  of  a  single  cell. 

3.  Survey  of  thermal  behavior  of  components  in 
lithium-ion  batteries 

Abuse  testing  of  batteries  gives  some  insight  into  failure 
mechanisms,  but  direct  characterization  of  the  materials  that 
constitute  the  battery  can  lead  to  an  understanding  of  what 
makes  a  battery  abuse  tolerant.  Below,  the  important 
exothermic  reactions  that  take  place  during  abuse  testing 
of  lithium-ion  batteries  are  summarized,  and  then  the  lit¬ 
erature  supporting  each  reaction  reviewed  in  order  to  extract 
quantitative  information  useful  for  modeling. 

The  following  exothermic  reactions  have  been  proposed 
to  take  place  during  abuse  testing  of  lithium-ion  batteries. 

(a)  Solid  electrolyte  interface  (SEI)  layer  decomposition: 
The  negative  electrode  is  protected  from  direct  reaction 
with  solvent  by  an  electronically  insulating  but 
ionically  conducting  film  called  the  “solid  electrolyte 
interface”.  This  layer  is  metastable  and  can  decompose 
exothermically  at  90-120  °C. 

(b)  Reaction  of  intercalated  lithium  with  electrolyte:  At 
elevated  temperatures  (>120  °C),  the  SEI  layer  does  not 
protect  the  negative  electrode  from  contact  with  the 
electrolyte  and  an  exothermic  reaction  between  inter¬ 
calated  lithium  and  electrolyte  can  occur. 

(c)  Reaction  of  intercalated  lithium  with  fluorinated 
binder:  Fluorinated  binder  can  react  exothermically 
with  lithiated  carbon. 

(d)  Electrolyte  decomposition:  The  electrolyte  can  decom¬ 
pose  exothermically  at  elevated  temperatures 
(>200  °C). 

(e)  Positive  active  material  decomposition:  In  the  oxidized 
state,  the  positive  active  material  can  decompose 
exothermically  and  give  off  oxygen.  This  oxygen  can 
react  exothermically  with  electrolyte.  Or,  perhaps  the 
positive  material  reacts  directly  with  electrolyte.  In  any 


case,  the  chemical  reduction  of  the  positive  active 
material  is  highly  exothermic. 

(f)  On  overcharge,  metallic  lithium  is  formed  which  can 
react  with  electrolyte. 

(g)  Lithium  metal  can  react  with  binder. 

(h)  Discharge  of  the  battery  releases  heat  due  to  entropy 
changes,  overpotentials,  and  ohmic  resistances. 

The  literature  supporting  each  of  these  exothermic  processes 
is  reviewed  below. 

3.1.  SEI  layer  decomposition 

Richard  and  Dahn  [7]  identifies  an  exothermic  peak  due 
to  SEI  decomposition  in  ARC  studies  of  carbon  anodes  in 
electrolyte  at  ^100  °C.  The  peak  is  relatively  insensitive  to 
the  amount  of  lithium  intercalated  into  the  carbon  and 
does  not  appear  if  fresh  carbon  is  heated  with  electrolyte. 
The  SEI  layer  is  assumed  to  consist  of  stable  (such  as 
Li2C03),  and  metastable  components  (such  as  (CH20- 
C02Li)2).  The  metastable  component  might  exothermi¬ 
cally  react  as  follows: 

(CH20C02Li)2  -►  Li2C03  +  C2H4  +  C02  +  \02  (la) 
or 

2Li  +  (CH20C02Li)2  -►  2Li2C03  +  C2H4  (lb) 

This  depletion  of  the  metastable  SEI  component  gives  rise  to 
a  peak  in  the  ARC  experiment.  The  composition  of  the  SEI 
depends  on  the  electrolyte  composition.  LiBF4-based  elec¬ 
trolytes  do  not  exhibit  a  peak  in  the  self-heating  rate  profile. 
With  LiPF6-based  electrolytes,  the  shape  and  behavior  of  the 
peak  depend  on  the  solvents  used.  For  example,  EC: diethyl 
carbonate  (DEC)  (1:1)  gives  a  shoulder  on  the  main  peak, 
whereas  EC: DEC  (1:2)  does  not. 

From  ARC  measurements  made  at  various  initial  tem¬ 
peratures  (plots  of  In  dT / d t  versus  1/7)  an  activation  energy 
of  1.35  kJ/mol  and  a  frequency  factor  of  1.7  x  1015  s1  were 
obtained.  The  ARC  peak  attributed  to  SEI  decomposition 
(1  M  LiPF6  in  EC:DEC  (1:2),  MCMB  carbon)  raises  the 
sample  temperature  about  15  °C,  which  corresponds  to 
257  J/g  of  carbon.  Increasing  the  surface  area  increased 
the  amount  of  SEI  and  correspondingly  the  magnitude  of 
the  associated  enthalpy. 

Dahn’s  group  has  examined  a  number  of  carbon  electrode 
materials  (fibers,  MCMB,  coke,  and  synthetic  graphites) 
using  ARC  [29]  and  found  the  initial  temperature  rise  is 
proportional  to  the  surface  area  of  the  carbon,  with  the 
exception  of  graphite. 

Maleki  et  al.  [3]  in  differential  scanning  calorimetry 
(DSC)  studies  of  lithiated  carbons  with  electrolyte  observed 
a  peak  at  ^100  °C  that  they  attributed  to  thermal  decom¬ 
position  of  the  SEI  layer. 

Zhang  et  al.  [2]  reported  DSC  results  for  mixtures  of 
lithiated  carbon  (MCMB-28)  and  electrolyte  (1  M  LiPF6  in 
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1:1  EC:DMC).  They  found  a  peak  beginning  at  —  130  °Cthat 
was  small  and  independent  of  the  degree  of  lithium  inter¬ 
calation;  the  peak  was  typically  —40  J/g,  though  a  maximum 
value  of  120  J/g  was  observed.  If  the  sample  was  70  wt.% 
active  material  and  8  wt.%  binder,  then  the  worst  case  heat 
of  reaction  is  186  J/g  of  carbon. 

3.2.  Reaction  of  intercalated  carbon  with  electrolyte 

Solvent  might  react  with  lithium  (either  metallic  or 
intercalated)  as  follows: 

2Li  +  C3H4O3  (EC)  -►  Li2C03  +  C2H4  (2a) 

2Li  +  C4H603  (PC)  ->  Li2C03  +  C3H6  (2b) 

2Li  +  C3H603  (DMC)  -►  Li2C03  +  C2H6  (2c) 

Roth  et  al.  [28]  reported  DSC  studies  of  lithiated  carbon 
anodes  (from  Sony  18,650  cells)  and  EC:PC:DEC/LiPF6 
electrolytes.  The  lithiated  carbon  anode,  when  dried, 
showed  peaks  above  300  °C.  Interestingly,  addition  of 
solvent  without  salt  had  no  effect  on  the  traces.  However, 
in  the  presence  of  electrolyte  (solvents  +  salt),  the  traces 
became  very  complex  with  peaks  at  —100,  —150,  —270, 
and  -300  °C. 

Richard  and  Dahn  [7]  studied  the  thermal  stability  of 
lithium  intercalated  graphite  in  electrolyte  by  differential 
scanning  calorimetry  and  ARC.  For  the  reaction  of  lithiated 
carbon  (MCMB)  and  electrolyte  (1  M  LiPF6  in  1:2 
EC:DEC),  they  report  an  enthalpy  of  1714  J/g  of  carbon. 
Since  Lio.gQ  was  used,  the  reaction  enthalpy  for  interca¬ 
lated  lithium  can  be  adjusted  to  1.54  x  105  J/mol  of  lithium. 
They  report  activation  energy  of  Ea  =  1.35  x  105  J/mol,  and 
a  frequency  factor  of  2.5  x  1013  s"1. 

Biensan  et  al.  [4]  report  thermal  studies  of  components 
from  lithium-ion  cells.  For  FiPF6  in  1:1:3  PC:EC:DMC 
electrolyte,  they  found  a  peak  at  —120  °C  (350  J/g)  attrib¬ 
uted  to  hixC6  +  electrolyte.  The  paper  does  not  specify  if 
the  heat  of  reaction  is  normalized  with  respect  to  the  carbon 
or  electrode  material,  or  electrode  material  and  electrolyte. 
Assuming  the  enthalpy  refers  to  grams  of  carbon,  and 
taking  the  capacity  of  the  carbon  as  330  mAh/g  gives 
2.84  x  104  J/mol  of  lithium.  The  use  of  FiBF4  and  FiTFSI 
salts  resulted  in  dramatically  different  DSC  traces.  How¬ 
ever,  the  enthalpies  for  the  reaction  with  lithiated  carbon 
were  not  reported. 

Fampe-Onnerud  et  al.  [11]  report  a  heat  of  reaction  of 
1497  J/g  for  the  reaction  of  FiC6  +  EC/DMC-FiPF6  with  an 
activation  energy  of  34  kJ/mol.  Assuming  a  capacity  of 
330  mAh/g,  this  gives  a  reaction  enthalpy  of  1.2  x  105  J/ 
mol  of  lithium. 

Du  Pasquier  et  al.  [18]  reported  the  enthalpies  for  the  sum 
of  SEI  decomposition  and  lithiated  carbon  (MCMB  25-28)/ 
electrolyte  (1  M  FiPF6  in  2:1  EC:DMC)  reactions.  The 
enthalpy  depended  on  the  degree  of  lithiation;  assuming  a 
capacity  of  330  mAh/g  for  MCMB  and  averaging  their 
values  gives  7.5±1.3xl04  J/mol  of  lithium. 


3.3.  Reaction  of  lithiated  carbon  with  fluorinated  binder 

Du  Pasquier  et  al.  [18]  suggested  that  lithium  (either 
metallic  or  intercalated)  can  react  with  a  fluorinated  binder 
as  follows: 

-CH2-CF2-  +  Li  -►  LiF  +  -CH=CF-  +\U2  (3) 

Alternatively,  fluorinated  binder  might  react  according  to 
this  scheme: 

2Li  +  RF2  —  2LiF  +  \  R2  (4) 

where  the  molecular  weight  of  RF2  is  64  g/(g  mol). 

Typically,  a  negative  electrode  contains  about  8-12  wt.% 
binder  with  the  balance  being  carbon.  If  the  electrode  is 
8  wt.%  binder,  the  carbon  has  a  capacity  of  330  mAh/g,  then 
a  fully  charged  electrode  contains 

(330  mAh/g  C6)  ( 0.92  g  C6/g  electrode) 

x  ( 3 .6  A  s/mAh)  (1/96, 500)  (mol  Li/ A  s) 

=  0.01 13  mol  Li/g  electrode 

(0.08  g  PVDF  /  g  electrode)  (2  mol  F  /  mol  PVDF) 
x  ( 1/64)  (mol  PVDF/g  PVDF) 

=  0.0025  mol  F/g  electrode 

So  a  fully-charged  negative  contains  at  least  four  times  as 
many  moles  of  Li  than  moles  of  fluorine.  At  first  sight,  this 
excess  of  lithium  seems  inconsistent  with  lithium  fluoride 
formation  as  the  reaction  product,  since  DSC  studies  [10] 
with  8  wt.%  PVDF  binder  indicate  the  heat  of  LixC6/binder 
reaction  increases  linearly  with  degree  of  lithiation  of  the 
negative.  On  the  other  hand,  other  workers  [4]  report  that  the 
heat  of  reaction  increases  linearly  with  the  amount  of  PVDF 
binder  from  10  to  15%  binder.  This  discrepancy  might  be 
due  to  the  availability  of  lithium  during  the  DSC  experiment. 
Since  lithium  must  diffuse  through  graphite  to  reach  the 
PVDF  binder,  not  all  the  lithium  is  available  for  reaction 
(indeed,  the  Li  and  the  PVDF  are  not  in  the  same  place, 
initially  at  least,  so  the  mere  stoichiometric  ratio  is  not  very 
meaningful).  As  the  degree  of  lithiation  increases,  more 
lithium  becomes  available  for  reaction.  On  the  other  hand, 
when  the  carbon  is  fully  lithiated,  the  reaction  with  PVDF  is 
limited  by  PVDF. 

Du  Pasquier  et  al.  [18]  reported  the  enthalpies  for  the 
reaction  of  lithiated  carbon  and  lithium  with  various  fluori¬ 
nated  binders.  The  authors  suggest  that  LiC6  reacts  at  lower 
temperatures  than  lithium  metal  due  to  the  higher  surface 
area  of  the  carbons.  They  found,  through  X-ray  analysis  of 
the  reaction  products,  that  LiF  is  the  major  inorganic  pro¬ 
duct.  For  PVDF  homopolymer,  assuming  a  hixC6  capacity  of 
330  mAh/g  gives  1.49  x  105  J/mol  of  LiC6. 

Maleki  et  al.  [3]  report  a  value  of  317  J/g  for  the  reaction 
of  hixC6  with  PVDF;  the  reaction  begins  at  —200  °C  and 
reaches  a  maximum  at  287  °C.  Later,  Maleki  et  al.  [10] 
reported  a  more  in-depth  study  of  the  lithiated  carbon/binder 
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reaction  in  which  various  carbons  and  various  binders  were 
examined.  The  authors  explain  that  hard  carbon  (HC-1227) 
produces  an  exothermic  DSC  peak  even  when  uncharged 
due  to  reaction  of  “irreversible  lithium”  with  PVDF.  The 
irreversible  lithium  is  produced  during  SEI  formation  on  the 
first  charging  cycle.  It  is  interesting  to  note  that  the  SFG-44 
produces  more  heat  than  the  MCMB  even  at  the  same  state 
of  charge.  Since  SFG-44  has  a  higher  surface  area  than 
MCMB,  this  is  consistent  with  the  proposal  that  the  reaction 
is  diffusion-controlled.  In  any  case,  using  the  slopes  from 
plots  of  heat  of  reaction  versus  charge,  gives  7.2  x  104  J/ 
mol  Li  (SFG-44),  5.9  x  104  J/mol  Li  (MCMB),  and  3.8  x 
104  J/mol  Li  (HC-1227). 

Biensan  et  al.  [4]  report  a  value  of  1500  J/g  for  the 
reaction  of  LixC6  with  PVDF;  the  reaction  begins  at 
~240  °C,  peaks  at  ~290  °C,  and  is  complete  at  ^350  °C. 

3.4.  Electrolyte  decomposition 

Published  work  on  electrolyte  decomposition  is  some¬ 
what  contradictory  and,  in  any  case,  the  situation  is  com¬ 
plicated  by  interactions  of  the  electrolyte  and  the  cathode. 
For  example,  ethylene  carbonate  decomposes  to  carbon 
dioxide  and  ethylene  oxide;  ethylene  oxide  can  exothermi¬ 
cally  decompose.  The  decomposition  of  EC  is  catalyzed  by 
lithium  nickelate  [13].  Other  workers  have  found  that  the 
amount  of  heat  released  in  the  reaction  between  Li0  5Co02 
and  electrolyte  is  not  affected  by  amount  of  electrolyte, 
within  certain  limits,  but  the  amount  of  heat  released  in  the 
reaction  between  Li^Mn204  and  electrolyte  increases  with 
increasing  electrolyte  [16]. 

Sloop  et  al.  [30]  studied  the  stability  of  LiPF6-EC:DMC 
electrolytes  at  85  °C.  Their  results  indicate  that  PF5  gas,  in 
equilibrium  with  LiPF6  (LiPF6  — >  LiF  +  PF5)  reacts  prefer¬ 
entially  with  EC  over  DMC  to  form  soluble  oligomers  of 
ethylene  carbonate  as  well  as  insoluble  precipitates  (perhaps 
phosphate  esters). 

Botte  et  al.  [9]  studied  the  thermal  stability  of  LiPF6- 
EC:EMC  electrolytes.  They  found  the  thermal  behavior  was 
very  sensitive  to  the  electrolyte  composition.  For  EC:EMC 
mixtures  they  report  DSC  traces  that  exhibit  a  slight 
endotherm  before  a  broad  exotherm.  For  1  M  LiPF6  solu¬ 
tions,  they  report  an  enthalpy  of  ~300  J/g  of  electrolyte.  The 
solvent  concentration  is  taken  as  0.0087  mol/g,  which  gives 
3.43  x  104  J/mol  of  solvent. 

Biensan  et  al.  [4]  report  a  heat  of  reaction  of  250  J/g  for 
decomposition  of  LiPF6-PC:EC:3DMC  electrolyte. 

Kawamura  et  al.  [43]  report  DSC  studies  of  various 
solvents.  They  found  that  diethyl  carbonate  was  more  reac¬ 
tive  than  dimethyl  carbonate  in  either  propylene  carbonate 
or  ethylene  carbonate  and  LiPF6  or  LiC104.  For  1  M  LiPF6 
electrolytes,  they  found  exotherms  between  230  and  280  °C 
with  heats  of  reaction  ~375  J/g  (1:1  PC:DMC  or  EC:DMC) 
and  515  J/g  (1:1  PC:DEC  or  EC:DEC).  Addition  of  1  wt.% 
water  to  the  electrolytes  reduced  the  magnitude  of  the 
exotherms  by  about  a  factor  of  two. 


3.5.  Positive  active  material  decomposition 

Charged  positive  active  materials  can  disproportionate  at 


elevated  temperatures  as  follows  [33]: 

Nio.8Coo.2CL  — >  ^Ni2.4Coo.604  +  ^02  (5a) 

Mn204  — >  Mn203  +  (5b) 

The  oxygen  released  might  react  with  solvent  as  follows: 
§02  +  C3H403  (EC)  -►  3C02  +  2H20  (6a) 

402  +  C4H603  (PC)  -►  4C02  +  3H20  (6b) 

302  +  C3H603  (DMC)  ->  3C02  +  3H20  (6c) 


The  work  of  Arai  et  al.  [44]  suggests,  for  delithiated  lithium 
nickelate,  that  oxygen  released  by  decomposition  of  \AX_ 
Ni02  completely  combusts  solvent  according  to  reaction  (6). 

MacNeil  and  Dahn  [8]  studied  the  reaction  kinetics  for 
LixCo02  in  non-aqueous  electrolyte.  They  report  a  good  fit 
with  their  data  using  a  reaction  enthalpy  of  265  J/g  of  LiCo02; 
multiplying  by  the  molecular  weight  of  LiCo02  and  assuming 
0.5  eq.  of  Co4+/mol  of  Li0.5CoO2  gives  5.18  x  104  J/eq.  They 
report  an  activation  energy  of  1.235  x  105  J/eq.,  and  a  fre¬ 
quency  factor  of  6.67  x  1011  s_1. 

Biensan  et  al.  [4]  report  enthalpies  for  several  different 
electrolyte/positive  electrode  reactions;  for  Li0.45CoO2, 
450  J/g  are  released  between  220  and  500  °C. 

Lampe-Onnerud  et  al.  [11]  report  the  heat  of  reaction  for 
Li0  5Co02  +  EC/DMC-LiPF6  occurs  with  two  DSC  peaks 
(total  615  J/g).  The  first  exothermic  process  has  an  activation 
energy  of  70  k J/mol. 

Venkatachalapathy  et  al.  [1]  reported  ARC/DSC  studies  in 
which  they  found  AH  of  641.7  J/g  for  LixNi0.8Co0.2O2  and 
381  J/g  for  LixCo02.  Unfortunately,  the  stoichiometry  (x)  of 
the  materials  was  not  disclosed.  These  authors  also  reported 
the  thermal  behavior  of  LiNi0.8Co0.2O2  in  18,650-size 
lithium-ion  cells  [6].  DSC  studies  of  the  individual  fully- 
charged  electrodes  taken  from  the  cells  show  two  exotherms 
(one  minor  and  one  large)  for  the  negative  and  one  major 
exotherm  for  the  positive  electrode.  The  AH  for  the  carbon 
electrode  is  317  and  641  J/g  for  the  LiNi0.8Co0.2O2  elec¬ 
trode.  Assuming  the  carbon  electrode  is  90%  active  material 
and  using  72  g/(g  mol)  for  carbon  gives  3.17  x  104  J/eq.  for 
the  LiC6/solvent  reaction.  Similarly,  assuming  the  LiNi0.8- 
Co0.2O2  electrode  contains  85%  active  material  and  196  g/ 
eq.  gives  1.48  x  105  J/eq.  for  the  LiNi0.8Co0.2O2/solvent 
reaction. 

Maleki  et  al.  [3]  reported  the  DSC  behavior  of  positive 
electrodes  containing  Sn-doped  LiCo02.  Washed  samples 
showed  an  exotherm  (146  J/g)  stretching  from  178  to 
250  °C.  Unwashed  samples  showed  an  exotherm  (407  J/g) 
beginning  at  167  °C  and  reaching  a  peak  at  214  °C. 

Lu  et  al.  [6]  reported  the  thermal  behavior  of  delithiated 
Lii  JCNi02.  They  found  the  AH  for  thermal  decomposition 
increased  linearly  with  increasing  v  from  v  =  0-0.5  to  a 
maximum  of  ~1 15  J/g;  beyond  v  =  0.5  the  enthalpy  actually 
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decreased  slightly  with  increasing  x.  For  v  =  0.5,  the  enthalpy 
corresponds  to  2.25  x  104  J/eq.  of  Ni02. 

Zhang  et  al.  [2]  carried  out  DSC  studies  of  positive 
electrode  materials  with  electrolyte  (1  M  LiPF6  in  1:2 
EC:DMC).  They  found  the  heat  generated  increased  almost 
linearly  with  degree  of  delithiation.  The  peak  temperature  of 
the  DSC  exotherms  shifted  to  lower  temperatures  as  the 
degree  of  delithiation  increased.  Their  results  indicate  that 
dried  samples  (electrolyte  removed)  released  significantly 
less  heat,  but  quantitative  results  were  not  provided. 

MacNeil  et  al.  [32]  studied  the  reaction  of  LixMn204  and 
electrolyte.  They  found  the  thermal  stability  of  LixMn204 
decreased  as  the  concentration  of  LiPF6  (in  EC/DEC) 
increased.  For  reaction  of  0.1  g  of  Mn204  with  0.1  g  of 
electrolyte,  they  found  a  heat  release  of  64  J  versus  a 
theoretical  (based  on  complete  oxidation  of  EC)  of  263  J; 
indicating  the  solvent  is  not  completely  oxidized. 

The  heat  of  reaction  for  positive  decomposition  is  coupled 
with  that  of  solvent  oxidation.  However,  the  mechanism  is 
not  clear.  Positive  active  material  could  react  directly  with 
solvent,  or  disproportionate  to  give  off  oxygen  that  reacts 
with  solvent. 

3.6.  Lithium  metal  reactions 

Little  data  is  available  on  these  reactions.  Work  by  von 
Sacken  et  al.  [15]  found  that  lithium  metal  anodes  were  less 
stable  than  lithiated  carbon  anodes.  There  is  some  evidence 
[26,28]  that  the  electrolyte  reaction  exhibits  a  peak  at  the 
melting  point  of  lithium  (~180  °C).  DSC  studies  with 
lithium  metal  and  PVDF  indicate  that  this  reaction  takes 
place  at  higher  temperatures  than  the  lithiated  carbon/PVDF 
reaction  [18,28].  However,  this  result  is  likely  due  to  the  low 
surface  area  of  the  materials.  In  a  battery,  lithium  metal  and 
PVDF  would  be  well  dispersed  and  possibly  in  intimate 
contact. 

Kawamura  et  al.  [43]  found  using  DSC  that  lithium  metal 
exhibits  an  exothermic  peak  with  dry  electrolytes  (1  M 
LiPF6  in  1:1  PC:DMC  or  1:1  EC:DMC  or  1:1  EC:DEC) 
at  ~180  °C.  They  suggest  that  the  SEI  layer  on  lithium 
protects  the  metal  from  reaction  until  the  melting  point. 
Addition  of  1  wt.%  water  decreased  the  temperature  of  the 
exothermic  peak  to  ~  140  °C  which  the  authors  attribute  to  a 
less  stable  SEI. 

3. 7.  Discharge  of  the  battery  releases  heat  due  to 
entropy  changes,  overpotentials,  and  ohmic  resistances 

The  amount  of  heat  released  due  to  entropy  changes  in 
lithium-ion  cells  is  small  relative  to  ohmic  heating.  For 
example,  Kim  et  al.  [31]  characterized  heat  generation 
of  LiJVln204  spinel  in  a  microcalorimeter.  At  low  rates 
(<1C)  of  charge/discharge,  the  thermal  behavior  reflected 
reversible  processes.  At  and  above  1C,  the  thermal  behavior 
was  dominated  by  irreversible  processes  (ohmic  and  over¬ 
potentials). 


3.8.  Summary  of  component  studies 

Table  la-d  summarizes  data  on  heat  release  from  thermal 
reactions  taking  place  in  lithium-ion  cells. 

The  model  requires  expressing  the  heats  of  reaction  in 
terms  of  moles.  Table  2  summarizes  the  heats  of  reactions  on 
a  molar  basis. 


4.  Estimation  of  reaction  parameters 

All  the  reactions  are  assumed  to  follow  an  Arrhenius 
expression  for  the  kinetic  rate  constant.  There  is  very 
limited  information  on  activation  energies  ( Ea )  and  fre¬ 
quency  factors  (k0).  In  order  to  estimate  these  values,  the 
following  procedure  was  followed.  First,  assuming  all  the 
reactions  are  pseudo- first-order,1  allowed  writing  a  rate 
expression  involving  activation  energy  and  frequency  fac¬ 
tor.  This  rate  expression  was  used  in  an  energy  balance  with 
two  unknowns:  the  activation  energy  and  frequency  factor. 
By,  somewhat  arbitrarily,  specifying  the  adiabatic  self¬ 
heating  rate  as  0.001  °C/min  at  the  onset  temperature  for 
thermal  runaway  (estimated  from  literature  review  or 
guessed,  see  Table  3),  the  activation  energy  and  frequency 
factor  could  be  computed.  The  results  are  summarized  in 
Table  3. 

These  kinetic  parameters  can  be  used  to  simulate  DSC  and 
ARC  experiments.  Figs.  1  and  2  show  the  simulation  results. 
The  DSC  curves  (Fig.  1)  indicate  the  magnitudes  of  the  heat 
released  per  gram  of  reaction  mixture  (1:1  molar  ratio  for 
binary  reactions).  The  reactions  involving  lithium  metal  are 
more  pronounced  than  the  analogous  LiC6  reactions  because 
of  the  low  MW  of  lithium.  The  DSC  traces  suggest  a 
sequence  for  the  reactions,  but  the  ARC  simulations  more 
clearly  illustrate  how  the  reactions  overlap. 

The  ARC  simulations  (Fig.  2)  indicate  the  SEI  decom¬ 
position  will  begin  first  at  70  °C.  At  85  °C,  the  lithium/ 
solvent  reaction  begins  and  then,  at  110  °C  the  negative/ 
solvent  reaction  begins.  The  LiC6/binder  reaction  begins  at 
160  °C  but  goes  at  a  lower  rate  than  the  NiCo02  decom¬ 
position  reaction,  so  by  170  °C  the  NiCo02  reaction  is 
taking  place  at  a  higher  rate  than  the  negative/binder  reac¬ 
tion.  The  solvent  decomposition  and  lithium/binder  reac¬ 
tions  both  begin  at  180  °C,  but  the  self-heating  rate  of 
solvent  decomposition  reaction  is  slightly  slower.  Finally, 
at  190  °C  decomposition  of  Mn204  begins.  Because  of  the 
overlapping  reactions,  the  role  of  individual  components  in 
thermal  runaway  is  not  clear;  this  is  exemplified  by  the 
reaction  between  binder  and  lithiated  carbon. 

The  PVDF  binder  reacts  only  with  the  negative,  if  there  is 
negative  reactant  left  after  reacting  with  solvent.  Even  then, 
the  positive  electrode  reaction  quickly  overwhelms  the 


1  The  assumption  of  first-order  reactions  is  only  made  in  estimation  of 
parameters.  In  the  battery  simulation  model  calculations,  second-order  and 
non-integral  order  reactions  are  considered  (see  Table  4). 
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Table  1 

Summary  of  reaction  parameters 


Reference 

SEI  Reaction 

LiC6/solvent 

LiC6/PVDF 

Energy  (J/g) 

Temperature 

Energy  (J/g) 

Temperature 

Energy  (J/g) 

Temperature 

range  (°C) 

range  (°C) 

range  (°C) 

(a)  Negative  electrode  (J/g  of  carbon,  fully  charged) 

[4] 

350 

110-150 

1500 

240-350 

[11] 

1497 

90-290 

[7,12,29] 

257 

100 

1714 

[3] 

100 

-1100 

220-300 

[2] 

186 

130 

[17] 

1460 

206-215 

[18]  (Li0.9C6) 

720 

110-260 

1394 

270-400 

(b)  Positive  electrode  decomposition  (per  g  of  active,  fully  charged) 

Reference 

LiCo02 

LiMn204 

LixNi0.8Co0.2O2 

Energy  (J/g) 

Temperature 

Energy  (J/g) 

Temperature 

Energy  (J/g) 

Temperature 

range  (°C) 

range  (°C) 

range  (°C) 

[5] 

115 

175-340b 

[3] 

146 

178-250 

(c)  Positive/solvent  reaction  (per  g  of  active) 

Reference 

LiCo02 

LiMn204 

LixNi0.8Co0.2O2 

Energy  (J/g) 

Temperature 

Energy  (J/g) 

Temperature 

Energy  (J/g) 

Temperature 

range  (°C) 

range  (°C) 

range  (°C) 

[4] 

450 

220-500 

450 

150-300 

600a 

180-500a 

[1,6] 

381 

235 

641.7 

216 

[11] 

615 

230-300 

[8] 

265 

170-240 

[2] 

625 

200-260 

350 

200-400 

1256 

200-230c 

[3] 

407 

167-250 

(d)  Electrolyte  (per  g  of  electrolyte) 

Reference 

Solvent 

Energy  (J/g) 

Temperature  range  (°C) 

[9] 

-285 

250-350 

[4] 

250 

130-220 

[11] 

155 

225-300 

[43] 

370-530 

230-280 

aLi0.3NiO2. 
b  Lio.5Ni02. 

c  Lio.25Ni02. 

Table  2 

Selected  heats  of  reactions  (per  mole) 

Reaction 

Equation 

Literature  values  (J/g) 

g/eq. 

J/eq.  or  J/mol 

SEI 

SEI— >  I 

186-257 

81.2a 

1. 5-2.1  x  104 

LiCVsolvent 

LiC6  +  dnsRS  — >  I 

350-1714 

81.2a 

0.28-1.39  x  105 

LiCVbinder 

LiC6  +  r>NBB  — >  I 

1100-1500 

81.2a 

0.89-1.2  x  105 

LiCo02 

265-625 

195.75 

0.52-1.2  x  105 

LixNi0.8Co0.2O2 

NiCo02  ->  I  +  \02 

600-1256 

139.8 

0.8-1. 8  x  105 

LiMn204 

Mn204  — >  I  +  402 

350-450 

180.8 

6.3-8. 1  x  104 

Solvent  decomposition 

S  — ►  I 

155-285 

100 

1. 5-2.9  x  104 

Solvent  oxidation 

02  T  DpsOS  — >  I 

1.43  x  106b 

Reactions  included  in  model.  I:  inert.  Mass  balance  kept  on  species  in  bold  type. 
a  Assuming  330  mAh/g. 

b  For  dimethylcarbonate  (from  CRC  Handbook,  63rd  ed.). 


R.  Spotnitz,  J.  Franklin/ Journal  of  Power  Sources  113  (2003)  81-100 


89 


Table  3 

Estimated  kinetic  parameters 


Reaction 

Tonset  (  L) 

Lpeak  (°C)  @  10  °C/min 

H/Cp  (K) 

Ea  (J/eq.) 

k0  (s  l) 

Li  +  ^libB  — >  I 

180 

280 

1.93  x  104 

2.86  x  105 

1.917  x  1025 

Li  +  »lisRS  — >  I 

85 

180 

1.02  x  104 

2.05  x  105 

9.41  x  1021 

LiC6  +  i>nbB  — >  C6  +  I 

160 

300 

3.71  x  103 

1.67  x  105 

1.79  x  1013 

LiCg  +  — >  C6  +  I 

110 

200 

9.81  x  102 

2.0  x  105 

1.95  x  1020 

Fi0.5NiCoO2  — >  I  +  ^02 

165 

225 

5.30  x  102 

3.94  x  105 

7.25  x  1039 

Mn204  — >  I  T  ^02 

190 

300 

7.23  x  102 

2.18  x  105 

1.06  x  1018 

S  ->  I 

180 

250 

1.91  x  102 

2.74  x  105 

5.14  x  1025 

SEI  -f  I 

70 

110 

1.43  x  102 

2.81  x  105 

7.88  x  1036 

H :  heat  of  reaction,  Cp :  heat  capacity. 


1000 

100 

o> 

§10 

c? 

1 

0.1 

60  100  140  180  220  260  300  340  380 

T,°C 

Fig.  1.  Simulated  DSC  experiments  for  cell  components. 


—  SEI  Decomp. 

-e-  Li/Solvent 

-K-  LiC6/Solvent 

NIC0O2  Decomp. 

-e-  Li/Binder 

-t- Solvent  Decomp. 

-a-  LiC6/Binder 

-e-  Mn2Q4  Decomp. 

10°C/min  scan  rate 


negative/binder  reaction,  so  it  is  not  clear  that  the  binder 
reaction  is  important.  If  the  heat  generated  from  the  SEI  and 
LiC6/solvent  reactions,  and  perhaps,  electrical  heating,  is 
sufficient  to  initiate  decomposition  of  the  positive,  then  the 
reactivity  of  the  binder  with  Li  may  be  a  moot  point. 


5.  Base  case  development 

Two  lithium-ion  chemistries  are  in  competition  for  use  in 
hybrid  electric  vehicles  and  42- V  batteries:  LiNizCo!_z02 
[22]  and  LiMn204  [21].  LiNizCo!_z02-cathodes  provide 


Li/Binder  Solvent 


Fig.  2.  Simulated  ARC  experiments. 
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Fig.  3.  Geometry  of  cell. 

better  life,  but  are  more  expensive  than  LiMn204-based 
cathodes.  To  model  the  abuse  behavior  of  these  chemistries, 
some  educated  guesses  must  be  made  as  to  the  internal 
design  of  the  corresponding  batteries.  Specifically,  the 
dimensions  of  the  electrodes,  collectors,  and  separator  must 
be  known,  as  well  as  the  compositions.  With  these  para¬ 
meters,  a  unit  cell  (see  Fig.  3),  consisting  of  A1  collector, 
positive,  separator,  negative,  Cu  collector,  can  be  specified. 
This  unit  cell  can  be  used  to  develop  model  parameters. 
Given  the  thickness  of  the  prismatic  cell,  and  the  unit  cell 
properties  (see  Table  7),  the  simulation  model  can  be 
specified. 

Given  the  quantities  and  heat  capacities  of  the  materials 
present  in  each  type  of  cell,  and  the  heats  of  reaction,  the 
maximum  temperature  rise  for  the  overall  battery  due  to 
each  reaction  can  be  estimated.  This  calculation  assumes  the 
heat  of  reaction  is  uniformly  distributed  throughout  all 
components  of  the  battery  and  so  does  not  account  for  local 
heating  effects  (“hotspots”).  Further,  discharge  and  charge 


conditions  will  be  different.  On  discharge  the  active  material 
is  limited  by  the  initial  state  of  charge,  but,  on  charge, 
lithium  metal  is  deposited  and  the  amount  of  positive  active 
is  increased  (up  to  a  certain  point).  Figs.  4  and  5  illustrate  the 
maximum  temperature  rise  possible  due  to  each  exothermic 
process  in  the  cell.  The  maximum  temperature  rise  for  each 
process  during  discharge  could  correspond  to  an  oven  test, 
short-circuit,  crush,  or  nail  test;  whereas,  the  maximum 
temperature  rise  for  each  process  during  charge  could 
correspond  to  an  overcharge  test. 

Both  chemistries  exhibit  similar  behavior.  The  tempera¬ 
ture  rise  due  to  SEI  decomposition  is  very  small  (<2  °C).  The 
negative/solvent  reaction  can  raise  the  cell  temperature 
~100  °C.  The  biggest  difference  between  the  two  chemis¬ 
tries  is  on  charge.  Since  the  manganese  oxide  cannot  be 
overcharged  (excess  charge  goes  to  solvent  oxidation),  while 
the  nickel  cobalt  oxide  electrode  can  be  further  delithiated, 
the  manganese  cell  generates  a  greater  temperature  rise  due 
to  solvent  oxidation.  Also  on  charge,  reactions  with  lithium 
metal  lead  to  significant  temperature  rises. 

On  discharge,  if  all  the  electrical  energy  went  to  cell 
heating,  the  cell  temperature  would  increase  well  above 
500  °C. 

The  solvent  oxidation  and  positive  active  decomposition 
reactions  should  be  considered  as  coupled  (oxygen  from  the 
positive  oxidizes  solvent).  Solvent  oxidation  causes  a  sig¬ 
nificant  temperature  rise  (>100  °C),  but  this  must  be  taken 
with  some  skepticism,  as  some  data  suggest  the  heat  released 
from  the  positive/solvent  reaction  is  much  less  than  that 
calculated  from  combustion  of  solvent  with  oxygen  released 
from  the  positive. 

The  temperature  rise  due  to  the  negative/binder  and 
negative/solvent  are  similar,  even  though  the  negative/binder 


Fig.  4.  Maximum  temperature  rise  due  to  exothermic  processes  during  discharge. 
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Fig.  5.  Maximum  temperature  rise  due  to  exothermic  processes  during  charge. 


reaction  is  much  more  exothermic  because  there  are  much 
fewer  moles  of  binder  than  solvent.  A  key  question  is  the 
relative  rates  of  the  negative/solvent  and  negative/binder 
reactions.  While  published  worked  has  shown  that  fluori- 
nated  binders  react  more  exothermically  with  lithiated  car¬ 
bon  than  non-fluorinated  binders  [10],  the  relative  rates  of 
the  negative/binder  and  negative/solvent  reactions  have  not 
been  clearly  addressed.  Biensan  et  al.  [4]  do  give  a  plot  of 
heat  release  versus  time  for  these  reactions,  but  an  isotherm 
of  700  °C  is  selected  and  the  time  scale  does  not  resolve  the 
two  processes. 


6.  Survey  of  thermal  modeling  of  lithium-ion  cells 


A  significant  amount  of  work  has  been  carried  out  to 
develop  mathematical  models  for  the  thermal  design  of 
lithium-ion  cells;  only  a  few  representative  papers  will  be 
described  here.  Chen  and  Evans  [35]  considered  the  problem 
of  heat  generation  in  a  prismatic  stack  of  lithium-ion  cells 
using  a  simple  energy  balance: 


dr  d2r  d2r  8 2 

~dl  =  kxdF  +  kyW  +  kz^  +  q 


(7) 


where  T is  the  temperature;  t  time;  p  the  cell  density;  Cp  the 
heat  capacity;  k  a  thermal  conductivity;  x,  y,  and  z  represent 
spatial  dimensions;  and  q  is  the  generation  rate.  The  heat 
generation  rate  was  taken  as  uniform  throughout  the  cell  and 
determined  from  experiments  with  actual  cells.  The  thermal 
conductivity  was  computed  by  averaging  reported  and  esti¬ 


mated  properties  of  the  cell  components.  The  calculations 
indicate  that,  even  for  large  cells  (20  cm  x  20  cm  x  20  cm), 
the  cell  temperature  would  not  reach  extreme  values,  nor 
would  the  temperature  gradients  be  large.  Abuse  behavior 
was  simulated  by  setting  the  internal  heat  generation  rate  at 
the  center  of  the  cell  to  a  high  value  for  a  short  time  and 
following  the  temperature  profile  with  time;  these  calcula¬ 
tions  indicated  that  thermal  runaway  could  occur  in  the  event 
of  a  localized  internal  short-circuit.  Kanari  et  al.  [36]  used  an 
approach  similar  to  Chen  and  Evans  [35],  but  considered  a 
cylindrical  cell  and  compared  their  results  to  experimental 
data  obtained  with  a  Sony  18,650  cell.  The  heat  generation 
rate  (« q )  was  computed  based  on  the  following  equation: 

q  =  -I(Eemi  -  V)  +  I  (r  )  (8) 

where  I  is  the  cell  current,  E  the  cell  voltage,  and  Eem f  the 
equilibrium  cell  voltage.  The  computed  temperatures  of  the 
can  at  both  ends  and  at  the  mid-point  of  the  tube  agreed  well 
with  experimental  data  during  discharge  and  charge  (within 
2  K).  Al  Hallaj  et  al.  [33]  developed  a  simplified  one¬ 
dimensional  model  that  gave  good  agreement  with  experi¬ 
mental  measurements  of  the  temperature  of  a  Sony  18,650 
cell.  They  used  the  model  to  estimate  the  thermal  behavior  of 
larger  (10-100  Ah)  cells  and  found  that  at  high  discharge 
rates  (1C)  the  internal  temperature  increased  by  ^60  °C. 
Sato  [37]  compared  experimental  and  model  predictions  for 
the  temperature  rise  during  discharge  of  an  80  Ah  Li-ion  cell 
and  obtained  good  agreement.  Most  recently,  Funahashi  et  al. 
[38]  used  a  model  similar  to  that  of  Kanari  et  al.  [36]  to 
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compare  with  actual  data  for  a  250  Wh  cell.  They  found 
excellent  agreement  between  measured  and  computed  sur¬ 
face  temperatures  at  a  1C  discharge  up  to  ~60%  DOD 
(~25  °C  temperature  rise);  at  higher  DOD  the  model 
over-predicted  the  temperature  rise  (model  estimated 
~36  °C  increase  versus  measured  ~30  °C).  The  discrepancy 
is  most  likely  because  the  HE  data  obtained  from  small  cell 
testing  does  not  accurately  represent  the  large  cell  HE 
behavior.  The  large  cell  heats  up  as  it  discharges  and  the 
small  cell  HE  data  does  not  reflect  this.  Since  polarization 
decreases  with  increasing  temperature  (as  both  mass  trans¬ 
port  and  reaction  kinetics  increase  with  temperature)  the 
large  cell  generates  less  heat  than  the  model  predicts.  By 
modeling  the  polarization  due  to  mass  transport  and  charge 
transfer,  a  more  accurate  estimate  of  the  thermal  behavior 
might  be  obtained.  Models  have  been  developed  that  include 
estimates  for  ohmic  and  polarization  voltage  losses  and  these 
models  have  been  demonstrated  to  give  good  predictions  of 
current-voltage  behavior  of  small  lithium-ion  cells  [39], 
however,  comparisons  with  thermal  behavior  have  not  yet 
been  reported. 

The  response  of  batteries  to  oven  heating  (an  abuse  test) 
has  been  predicted  using  a  model  that  accounts  for  the 
kinetic  behavior  of  the  positive  and  negative  electrodes 
[37].  Hatchard  et  al.  [12]  obtained  reasonable  agreement 
between  model  calculations  and  experimental  oven  tests  for 
18,650-size  cells.  They  used  the  model  to  predict  the  effect 
of  increasing  the  radius  of  a  cylindrical  cell  and  of  increasing 
the  thickness  of  a  prismatic  cell.  In  both  cases,  they  found  a 
critical  size  at  which  thermal  runaway  would  result. 


7.  Model  development 

For  simplicity,  a  cell  is  modeled  as  a  slab  consisting  of 
layers  of  various  materials  (see  Fig.  3).  Only  one  dimension 
is  considered  (v-direction),  the  slab  is  considered  uniform  in 
the  y  and  z  directions.  By  symmetry,  only  half  of  the  cell  is 
considered.  The  origin  of  the  coordinate  system  is  at  the 
center  of  the  cell.  Using  a  control  volume  approach  [40], 
mass  and  energy  balances  are  carried  out  for  each  layer  of 
the  cell.  The  positive,  separator,  and  negative  layers  are 
composites.  The  electrodes  consist  of  active  material,  elec¬ 
trolyte,  binder  and  inert  materials,  while  the  separator  con¬ 
sists  of  electrolyte  and  inert  material;  shutdown  (melting)  of 
the  separator  is  not  considered  in  the  model.  The  boundary 
condition  at  the  centerline  is  that  the  heat  flux  is  zero,  and  the 
boundary  condition  at  the  can,  is  that  heat  can  transfer  by 
convection  and  radiation.  Table  4  summarizes  the  kinetic 
equations,  while  Table  5  lists  the  mass  balance  equations  and 
Table  6  lists  the  source  terms  in  the  energy  balance  (see 
Eq.  (7)).  The  set  of  differential  equations  was  solved  using 
the  IDA  routine  of  Hindmarsh  and  Taylor  [41]. 

The  forms  of  the  kinetic  expressions  in  Table  4  are 
explained  as  follows.  The  rates  of  the  positive  (7?P),  solvent 
( Rs )  and  SEI  (7?SEI)  decomposition  reactions  are  all  taken  as 


Table  4 

Kinetic  equations 

Reaction 

Equation 

Positive  decomposition 

Rp  =  kpnipo 

Solvent  decomposition 

Rs  =  kscs 

Negative  solvent 

Rns  =  Ln^ns,  rNS  =  ~^xsaSjNlNwNrS£N, 

Fe 

+  - 

mNaSyNpE 

Negative  binder 

RnB  =  &NB<2s,NWn/nX^n(1  —  £N)c5jWb 

Li  diffusion  in  negative 

Nn  =  —aSiNlN^Nkm,NC^(xs  —  Xb) 

SEI  decomposition 

RsEI  =  ^SEI^E 

Lithium  solvent 

D  ,  7  ^Li^Li 

^LiS  =  kusLN^N  1  .  ,  cs 

1  +  kumu 

Lithium  binder 

_  kumu  wb/n 

LlB  LlB  1  +  ^Li^Li  #MW 

Table  5 

Mass  balance  equations 

Species 

Symbol 

Mass  balance 

Positive  active  in 

mPO 

oxidized  state 

Positive  active 

mPR 

in  reduced  state 

Solvent  in  the  positive, 

Cs 

negative,  and  separator 
Oxidizable  solvent 

Co 

in  the  positive 

Reducible  solvent 

Cr 

in  the  negative 
Intercalated  lithium 

xb 

in  negative 

SEI  material 

mE 

dmpo 
d  t 

dmpR 

d  t 


i 

~F 


dcs 
d  t 


=  - Rs 


dco  _  —  v?sRp  +  vs(ioc/F) 
d  t  GpLp 

dcR  _  VnsRnS  +  ^LiS^LiS 
£n^N 

Nn 

Ln(1  —  £n)cnc^ 
RsEI 


d  t 
dxb 

d  mE 
d  t 


Li  metal 


dmu  p  p  iu 
mu  —7—  =  —Rus  ~  Rub  +  — 

d  t  F 


Binder  in  negative 


wB 


dwe 
d  t 


Bmw 

In 


(^LiB  +  Rnb) 


Table  6 

Energy  balance  source  terms 


Test  or  element 

Equation 

Oven 

N)hm  =  0 

Short 

^ohm  —  1  "  1  i  —  Fk^Xfo 

Oi 

Overcharge 

N)hm  =  inr 

Oi 

Nail 

Crush 

c  -2  H  |  -2  E  •  'Yh  i 

‘Johm  —  *  +  *nail  _  ’  Znail  —  AT 

v  i  &  nail  ’  shorted 

c  _  a  E  1  .-2  E  .  Y  i 

°ohm  —  1  ~r  f crush  _  j  fcrush  —  AJ 

C^cell  ^crushed 

Collector 

8  —  ^ohm  • 

Positive  active 

Sp  =  /Epos  +  RpHp  +  LposRsHs  +  —  Hp  +  Snhm 

Separator 

8s  =  i  Rsep  T  FsepRsHs  T  ^ohm 

Negative  active 

Sn  =  iEneg  +  R^sHns  +  RsHs  +  RnbHnb 

+RsElHsEI  +  ^LiS^LiS  +  ^LiB^LiB  +  N)hm 
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Table  7 

Physical  properties  of  base  case  unit  cells  [20,42] 


Cbase 

47base 

LiNiCo02 

Graphite 

LiMn204 

Graphite 

Weight  fraction 

Active 

0.84 

0.92 

0.87 

0.94 

Binder 

0.06 

0.08 

0.06 

0.06 

Carbon  black 

0.1 

0 

0.07 

0 

Solids  density  (g/cm3) 

4.584 

2.211 

4.202 

2.223 

Heat  capacity  (J/(g  K)) 

0.651 

0.720 

0.672 

0.641 

Thermal  conductivity 

0.061 

0.063 

0.062 

0.065 

(W/(cm  K)) 

Apparent  density  (g/cm3) 

2.75 

1.08 

2.62 

1.08 

Porosity 

0.400 

0.512 

0.377 

0.514 

Electrolyte-filled 

Density  (g/cm3) 

3.233 

1.697 

3.074 

1.701 

Heat  capacity  (J/(g  K)) 

0.823 

1.113 

0.839 

1.064 

Thermal  conductivity 

0.037 

0.032 

0.039 

0.033 

(W/(cm  K)) 

Weight  fraction 

Electrolyte 

0.149 

0.364 

0.148 

0.365 

Solids 

0.851 

0.636 

0.852 

0.635 

first-order;  this  seems  reasonable  given  the  limited  informa¬ 
tion  available  about  them.  The  rate  of  the  negative  solvent 
reaction  (RNS)  is  taken  to  be  inversely  proportional  to  the 
amount  of  SEI  material  and  proportional  to  the  surface  area 
of  the  negative  and  the  surface  concentration  of  lithiated 
carbon;  this  type  of  expression  is  similar  to  that  used  by 
Richard  and  Dahn  [7] .  The  surface  concentration  of  lithiated 
carbon  was  determined  using  a  lumped  mass  balance  and  the 
flux  expression  given  in  Table  4  for  lithium-ion  diffusion  in 
carbon.  The  reaction  of  lithiated  carbon  and  binder  was 
taken  to  be  proportional  to  the  amount  of  binder  and  the 
surface  concentration  of  lithium.  The  reactions  with  metallic 
lithium  were  taken  as  first-order  in  solvent  or  binder,  and  to 
depend  only  slightly  on  the  amount  of  lithium.  The  ku 
parameter  was  set  to  a  large  value,  to  effectively  make 
the  term  kumu/ (1  +  kuwiu)  equal  to  unity  when  mu  ^  0, 
and  zero  when  mu  =  0. 

The  physical  properties  of  the  materials  commonly  found 
in  lithium-ion  cells  were  taken  from  the  literature  [19]; 
Table  7  lists  the  physical  properties  used  and  Table  8  details 
the  design  parameters  for  the  base  case  cells. 

8.  Case  studies 

8.1.  Oven  test 

For  a  175  °C  oven  test,  Fig.  6  shows  how  the  temperature 
profile  develops  in  the  battery  over  time.  The  figure  shows 
the  temperature  distribution  in  one  half  of  the  battery  (the 
other  half  is  the  same  by  symmetry);  the  core  temperature  is 
shown  on  the  left  and  the  can  temperature  on  the  right.  As  the 
battery  heats  up,  the  temperature  is  highest  at  the  can  surface 


Table  8 

Design  parameters  for  base  case  unit  cells 


Parameter 

NiCo02/C  cell 

Mn2Q4/C  cell 

NiCo02 

Carbon 

Mn204 

Carbon 

Active 

Capacity  (mAh/g) 

160 

330 

115 

330 

Capacity  (mAh/cm2) 

1.4784 

1.4784 

1.2006 

1.2006 

Loading  (mg/cm2) 

11 

4.87 

12 

3.87 

Coat  thickness  (pm) 

40.0 

45.1 

45.8 

35.8 

Separator  thickness  (pm) 

25.4 

25.4 

40.0 

40.0 

Separator  porosity  (%) 

40.0 

40.0 

40.0 

40.0 

Collector  thicknes/2  (pm) 

10.0 

5.0 

10.0 

5.0 

Unit  cell  thickness  (pm) 

125.5 

136.6 

Total  area  (cm2) 

7730.4 

4997.5 

Unit  cell  area  (cm2) 

23.78 

23.78 

20.20 

20.20 

Overall  thickness  (cm) 

4.08 

3.38 

Overall  volume  (cm3) 

97.01 

68.29 

Number  of  unit  cells 

325 

325 

247 

247 

Coating  weight  (g) 

85.0 

37.6 

60.0 

19.3 

Coating  volume  (cm3) 

30.9 

34.9 

22.9 

17.9 

Active  weight  (g) 

71.43 

34.63 

52.17 

18.18 

Binder  weight  (g) 

5.10 

3.01 

3.60 

1.16 

Carbon  weight  (g) 

8.50 

0.00 

4.20 

0.00 

Electrolyte  weight  (g) 

14.93 

21.52 

10.40 

11.12 

Collector  weight  (g) 

20.78 

34.32 

13.43 

22.19 

Electrode  weight  (g) 

120.75 

93.49 

83.81 

52.65 

Separator  polymer  (g) 

10.84 

11.03 

Separator  electrolyte  (g) 

9.42 

9.60 

Wet  separator  weight  (g) 

20.26 

20.63 

Can  weight  (g) 

3.77 

3.20 

Jellyroll  weight  (g) 

234.50 

157.08 

and  steadily  decreases  toward  the  core  (see  4999.69  s).  As 
the  exothermic  reactions  are  activated  and  start  to  release 
heat,  the  interior  temperature  of  the  battery  increases,  and 
the  temperature  reaches  a  maximum  at  the  core  and  steadily 
decreases  toward  the  wall.  At  5839  s,  thermal  runaway 
ensues  starting  at  the  innermost  cell  of  the  battery.  The 
contribution  of  the  various  reactions  can  be  visualized  with 
the  aid  of  Fig.  7.  This  figure  shows  the  cumulative  energy 
released  from  each  reaction  as  a  function  of  time.  The  cause 


Fig.  6.  Temperature  profiles  for  oven  test  (175  °C)  with  NiCoQ2  case. 
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of  the  runaway  can  be  associated  with  the  negative/solvent 
reaction  (Q-RNS),  the  positive  decomposition  reaction  (Q- 
RP),  the  solvent  decomposition  reaction  (Q-RS),  and  pos¬ 
sibly  the  negative/binder  reaction  (Q-NB). 

Since  the  value  for  the  negative  surface  area  used  to  obtain 
Figs.  6  and  7  is  somewhat  low  (0.8  m2/g),  cases  were  run 
with  higher  values  (2.4  and  4.8  m 2/g).  The  temperature 
profiles  at  the  start  of  thermal  runaway  are  shown  in 
Fig.  8.  As  the  surface  area  of  the  negative  active  material 
increases,  the  time  for  thermal  runaway  to  initiate  decreases, 
and  the  position  in  the  cell  where  thermal  runaway  initiates 
moves  closer  to  the  wall.  These  results  can  be  explained  as 
follows.  When  the  negative  area  is  low  (0.8  m2/g),  the 
negative/solvent  and  negative/binder  reactions  rates  initially 
generate  heat  at  a  rate  low  enough  that  it  can  be  dissipated 
through  the  can  wall.  This  sets  up  a  temperature  profile  with 
the  maximum  temperature  at  the  core.  As  the  temperature 
increases,  the  rate  of  heat  generation  exceeds  the  rate  of  heat 
dissipation  and  runaway  ensues  starting  at  the  core  (highest 
temperature).  Increasing  the  negative  area  increases  the  rate 


Element  position  (from  core) 

Fig.  8.  Temperature  profiles  for  LiNiCo02  with  varying  negative  active 
areas  (175  °C  oven). 


of  heat  release  from  the  negative/solvent  and  negative/binder 
reactions.  In  this  case,  just  before  thermal  runaway,  the 
temperature  profile  is  such  that  the  highest  temperature  is 
at  the  can  wall.  So  the  heat  generation  rate  is  highest  at  the 
can  wall  and  lowest  at  the  core.  Since  the  heat  generated 
nearest  the  can  wall  is  dissipated  by  convection,  the  tem¬ 
perature  near  the  can  wall  becomes  lower  than  the  interior. 
The  interior  temperature  rises  due  to  heat  generation  and  the 
temperature  profile  shows  a  maximum  because  the  rate  of 
heat  generation  is  decreasing  toward  the  core  (lowest  tem¬ 
perature).  In  all  cases,  the  contribution  of  the  negative/binder 
reaction  appears  small  relative  to  that  of  the  negative/solvent 
reaction. 

To  see  if  the  negative/binder  reaction  contributed  to 
thermal  runaway,  a  case  was  run  using  a  negative  area  of 
4.8  m2/g  and  setting  the  enthalpy  for  the  negative/binder 
reaction  to  zero.  The  results  were  almost  identical  to  the  case 
where  the  enthalpy  of  the  negative/binder  reaction  was 
considered  (without  the  binder  reaction,  the  time  to  runaway 
increased  by  only  4.5  s).  Thus,  even  though  some  energy  is 
released  by  the  negative/binder  reaction,  it  is  insignificant 
compared  the  positive  and  solvent  decomposition  reactions. 
Under  the  conditions  described  here,  the  binder  does  not 
play  a  role  in  thermal  runaway. 

The  Mn204  cell  differs  from  the  NiCo02  cell  in  several 
important  ways.  Mn204  is  more  thermally  stable  than 
NiCo02.  The  Mn204  cell  is  thinner,  so  better  able  to 
dissipate  heat.  The  Mn204  cell  contains  less  binder  in  the 
negative  (6  wt.%  versus  8  wt.%  for  the  NiCo02  cell).  The 
difference  in  the  cells  reflects  differences  observed  in  com¬ 
mercial  cell  designs.  The  results  of  a  175  °C  oven  test  for  the 
Mn204  cell  are  shown  in  Figs.  9  and  10.  Thermal  runaway 
occurs  in  the  following  manner  (see  Fig.  10).  The  SEI 
reaction  contributes  a  small  amount  of  heat,  but  it  over¬ 
shadowed  by  the  negative/solvent  reaction.  The  solvent 
decomposition  reaction,  once  activated,  proceeds  to  com¬ 
pletion.  The  negative/binder  reaction  is  inhibited  by  the  slow 
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rate  of  lithium  diffusion  out  of  the  negative  electrode  and 
does  not  release  a  significant  amount  of  energy.  Thermal 
runaway  results  from  the  positive  electrode  decomposition 
reaction. 

In  the  Mn204  simulations,  the  solvent  decomposition 
reaction  contributes  a  significant  amount  of  energy.  If  a 
lower  porosity  negative  electrode  is  used  (~25%  versus  50% 
for  base  case),  the  amount  of  solvent  will  be  significantly 
reduced.  Figs.  11  and  12  show  the  results  of  reducing  the 
negative  porosity  and  using  a  higher  value  for  the  mass- 
transfer  coefficient  of  lithium  diffusion  in  the  negative 
(10-6  cm/s  versus  10-7  cm/s  for  base  case).  In  this  case 
the  battery  reaches  a  maximum  temperature  of  only 
~242  °C.  The  SEI  decomposition,  solvent  decomposition, 


negative/binder,  and  negative/solvent  reactions  proceed  to 
completion  without  raising  the  cell  temperature  to  the  point 
where  the  positive  electrode  reaction  runs  away.  However, 
thermal  runaway  does  result  if  the  surface  area  of  the 
negative  electrode  is  increased  to  4.8  m2/g,  because  the  cell 
temperature  rises  above  265  °C,  and  the  positive  electrode 
decomposition  reaction  takes  off. 

These  simulations  indicate  the  negative/binder  reaction  is 
relatively  unimportant  in  the  oven  test.  The  binder  must 
compete  with  the  more  facile  solvent  reaction  to  react  with 
the  negative  electrode.  Based  on  the  understanding  gener¬ 
ated  by  these  simulations,  it  is  possible  to  predict  that  the 
binder  reaction  could  be  important  in  situations  where  the 
amount  of  solvent  is  limited  (low  porosity  electrodes)  and  a 


Fig.  10.  Cumulative  energy  generated  vs.  time  for  Mn2Q4  cell  (175  °C  oven  test).  Negative  area  4.8  m2/g. 
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Fig.  11.  Temperature  versus  time  for  Mn204  cell  (175  °C  oven).  Negative  porosity  ~25%,  kmn  =  10  6  cm/s. 


high  binder  loading  (10  wt.%)  is  used.  However,  both  of 
these  conditions  are  contrary  to  producing  high-rate  negative 
electrodes. 

8.2.  Short-circuit  test 

In  this  test,  the  possibility  exists  for  electrical  heating  to 
raise  the  temperature  of  the  cell  to  the  point  where  thermal 
runaway  initiates.  The  current  is  limited  by  diffusion  of 
lithium  in  the  negative,  and  so  either  increasing  the  mass- 
transfer  coefficient  for  lithium  in  the  negative  or  increasing 
the  surface  area  of  the  negative  allows  for  higher  currents 
and  so  more  rapid  heating  of  the  battery.  Despite  the  rapid 
heating,  the  positive  decomposition  reaction  is  still  respon¬ 
sible  for  thermal  runaway.  Increasing  the  mass-transfer 
coefficient  with  the  negative  area  set  to  4.8  m2/g  reduces 
the  time  to  thermal  runaway  to  ~  13.5  s,  but  does  not  change 
the  mechanism  for  thermal  runaway. 


As  in  the  oven  test,  thermal  runaway  in  the  short-circuit 
test  is  a  result  of  activating  the  positive  decomposition 
reaction  by  heating  the  battery  to  a  high  temperature. 
However,  in  the  short-circuit  test,  the  battery  can  be  heated 
much  more  rapidly  by  irreversible  heat  generated  from 
passing  current  through  the  electrodes.  This  rapid  heating 
sets  up  a  steep  temperature  profile  with  the  highest  tem¬ 
perature  at  the  core,  where  the  positive  decomposition 
reaction  can  be  activated  before  the  other  exothermic  reac¬ 
tions  (such  as  the  solvent  reactions)  have  progressed  sig¬ 
nificantly.  This  contrasts  with  the  oven  test,  in  which  the 
battery  heats  up  much  more  slowly,  allowing  the  solvent 
reactions  to  progress  significantly  before  the  positive  elec¬ 
trode  is  activated.  Thus,  with  the  NiCo02  cell,  the  binder 
plays  an  even  less  important  role  in  thermal  runaway  during 
the  short-circuit  test  than  the  oven  test. 

With  the  manganese  cell,  the  positive  decomposition 
reaction  takes  place  at  a  higher  temperature  than  in  the 


Fig.  12.  Cumulative  energy  versus  time  for  Mn204  cell  (175  °C  oven).  Negative  porosity  ~25%,  kmn  =  10  6  cm/s. 
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Fig.  13.  Cumulative  energy  versus  time  for  short-circuit  test  with  Mn204  cell.  Mass-transfer  coefficient  for  Li  in  negative  10  7  cm/s  and  negative  area 
4.8  m2/g. 


position  (from  core) 


Fig.  14.  Temperature  profiles  for  NiCo02  overcharge  test. 


NiCo02  cell,  so  the  solvent  reactions  have  more  time  to 
proceed  (see  Fig.  13).  Discharge  of  the  negative  and  reaction 
of  lithiated  carbon  with  solvent  deplete  the  negative,  so  the 
binder  reaction  takes  place  to  only  a  very  limited  extent. 
Still,  the  cell  goes  into  thermal  runaway  due  to  activation  of 
the  positive  decomposition  reaction. 

As  in  the  oven  test,  the  binder  does  not  compete  effec¬ 
tively  with  solvent  to  react  with  lithiated  carbon,  and  so 
plays  a  limited  role  in  the  thermal  runaway  behavior. 

8.3.  Overcharge  test 

Overcharging  a  lithium-ion  battery  can  lead  to  lithium 
metal  formation  which  opens  the  possibilities  for  further 
reactions  with  the  solvent  and  binder. 

The  temperature  profiles  for  an  overcharge  test  (1  mA / 
cm2,  ~C!  1.7  rate)  with  the  base  case  NiCo02  battery  are 
shown  in  Fig.  14  and  the  cumulative  energy  source  terms  in 
Fig.  15.  The  cell  heats  up  slowly  for  the  first  ~5400  s,  at 


which  point  the  electrodes  are  fully  charged  and  the  reac¬ 
tions  shift  from  intercalation  at  the  negative  to  deposition  of 
lithium  metal,  and  from  lithium  de-intercalation  at  the 
positive  to  oxidation  of  solvent.  The  solvent  oxidation 
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Fig.  16.  Temperature  profiles  for  Mn204  overcharge  test.  Negative  area 
4.8  m2/g. 

reaction  due  to  overcharge  (Q-P-OC  term  in  Fig.  15)  gen¬ 
erates  heat  at  a  much  higher  rate  than  the  irreversible  heating 
due  to  the  intercalation  reactions.  As  the  temperature  rises, 
first  the  lithium/solvent  reaction,  then  the  negative/solvent 
reaction  become  activated,  and  thermal  runaway  occurs. 
Unfortunately,  the  simulation  fails  at  this  point  because 
the  gradients  become  too  steep  to  follow.  When  thermal 
runaway  occurs,  neither  the  lithium/binder  reaction  nor  the 
positive  decomposition  reaction  are  contributing  significant 
energy.  Runaway  is  mainly  due  to  metallic  lithium  reacting 
with  solvent. 

An  overcharge  test  with  the  base  case  Mn204  cell  (1  mA/ 
cm2,  r^C/1.6  rate)  yields  similar  results  as  with  the  NiCo02 
battery  (see  Figs.  16  and  17).  At  about  900  s,  the  positive 
electrode  is  fully  charged  and  the  reaction  shifts  to  solvent 
oxidation.  The  solvent  oxidation  reaction  heats  the  battery 
and  activates  the  lithium/solvent  reaction.  This  further 
increases  the  rate  of  battery  heating  and  activates  the  nega¬ 
tive/solvent  reaction.  With  thermal  runaway  well  in  pro¬ 
gress,  the  solvent  decomposition  and  lithium/binder 
reactions  are  activated.  To  verify  the  binder  plays  no  role, 
the  simulation  was  repeated  with  the  enthalpies  of  the  binder 


Fig.  17.  Heat  source  terms  for  Mn204  overcharge  test.  Negative  area 
4.8  m2/g. 


reactions  set  to  zero,  and  virtually  identical  results  were 
obtained. 

The  major  experimental  observation  with  the  overcharge 
test  is  that  runaway  is  sensitive  to  the  charging  current.  The 
simulation  model  shows  the  same  behavior.  A  simulation 
run  (Mn204  cell)  with  a  charging  current  of  0.5  mA/cm2 
does  not  result  in  thermal  runaway. 

If  the  solvent  is  depleted,  then  the  binder  might  play  a 
role.  The  stoichiometric  numbers  for  the  electrolyte/lithium 
and  electrolyte/binder  reactions  were  changed  from  0.5  to 
2.5  to  exaggerate  the  effect.  However,  the  battery  did  not  get 
hot  enough  to  activate  the  binder  reactions. 

8.4.  Nail  test 

This  model,  being  one-dimensional,  can  only  very  crudely 
approximate  the  nail  test.  In  this  model,  the  heat  released 
from  the  nail  is  instantaneously  distributed  in  the  plane 
perpendicular  to  the  nail  movement.  In  an  actual  nail  test, 
a  local  hotspot  could  form  at  a  particular  place  in  the  plane  of 
the  battery.  A  hotspot  will  form  because  the  heat  generated 
in  the  nail  cannot  be  dissipated  as  fast  as  it  is  generated  even 
though  the  thermal  conductivity  of  batteries  is  typically  10  x 
higher  in  the  plane  of  the  electrodes  (due  to  the  metal  current 
collectors)  than  in  the  direction  the  nail  penetrates  [20] .  To 
simulate  a  hotspot,  since  the  model  does  not  account  for 
thermal  gradients  in  the  plane  of  the  electrodes,  the  heat 
generation  in  the  nail  was  set  to  a  large  value.  If  the  nail 
moves  at  8  cm/s,  all  the  cells  are  rapidly  shorted  and  the 
results  are  similar  to  a  short-circuit  test;  thermal  runaway 
starts  at  the  innermost  cell  and  is  due  to  the  positive 
decomposition  reaction.  If  a  much  lower  nail  speed 
(0.01  cm/s)  is  used,  thermal  runaway  initiates  at  the  outer¬ 
most  cell;  the  outermost  cell  heats  up  from  the  current 
flowing  through  the  nail,  which  activates  the  positive  decom¬ 
position  reaction. 

8.5.  Crush  test 

In  practice,  this  test  is  carried  out  in  a  variety  of  ways.  The 
cell  can  be  crushed  with  a  round  bar  until  the  cell  voltage 
drops  to  a  low  value,  or  the  cell  can  be  crushed  by  a  certain 
fraction  of  its  thickness  (for  example,  to  50%  of  its  original 
thickness).  To  simulate  the  effect  of  crush,  one  or  more  cells 
is  shorted  so  that  the  current  generated  by  the  battery  flows 
through  those  cells.  As  with  the  simulated  nail  test,  this 
simulation  suffers  by  using  a  one-dimensional  model.  In 
reality,  the  short  formed  from  a  crush  would  be  localized  to  a 
spot  on  the  electrode  while  the  one-dimensional  model  must 
assume  the  short  is  uniformly  distributed  across  the  plane  of 
the  electrode.  Since  Joule  heating  is  inversely  proportional 
to  the  area  traversed  by  the  current,  the  one-dimensional 
model  vastly  underestimates  the  amount  of  heat  generated. 

Setting  the  model  so  that  the  innermost  cell  is  shorted 
gives  results  almost  identical  to  those  of  a  short-circuit  test. 
The  battery  heats  up  due  to  overpotential  losses  at  the 
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electrodes  and  activates  the  negative/solvent  reaction.  This 
further  heats  the  battery,  activating  the  positive  decomposi¬ 
tion  reaction  and  causing  thermal  runaway. 

9.  Conclusions  and  recommendations 

The  simulations  indicate  that  the  negative  electrode  bin¬ 
der  plays  a  relatively  unimportant  role  in  thermal  runaway. 
The  amount  of  binder  is  limited,  and  the  binder  reaction 
must  compete  with  the  more  facile  solvent  reaction  for 
lithiated  carbon.  Reactions  involving  lithiated  carbon  are 
somewhat  hindered  because  the  lithium  must  diffuse 
through  the  solid  carbon  to  react  at  the  surface.  When  the 
battery  is  overcharged,  substantial  amounts  of  lithium  metal 
can  form.  The  lithium  may  be  in  intimate  contact  with  the 
binder  and  can  react  directly  without  being  hindered  by 
solid-state  diffusion.  Still,  the  solvent  reaction  is  preferred 
and  can  induce  thermal  runaway  without  the  aid  of  the 
lithium/binder  reaction. 

The  short-circuit  test,  and  especially  the  nail  and  crush 
tests,  can  cause  rapid  and  localized  heating  of  the  battery. 
This  localized  heating  results  in  high  temperatures  that 
activate  all  the  chemical  reactions  in  the  cell.  However, 
since  the  relative  rates  of  the  reactions  do  not  change,  the 
binder  reaction  still  plays  a  minor  role. 

The  simulations  reported  here  were  unable  to  rationalize 
the  experimental  observation  that  fluorinated  binders  fail  in 
the  nail  test  above  4.0  V,  whereas  non-fluorinated  binders 
pass  the  nail  test  up  to  4.3  V  [4].  One  possibility  is  that  at 
very  high  rates,  the  binder/negative  reaction  overtakes  that 
of  the  electrolyte/negative  reaction.  Also,  the  simulations  do 
not  account  for  the  possibility  that  the  battery  can  eject 
electrolyte.  If  the  electrolyte  is  displaced,  then  it  would  no 
longer  compete  with  the  binder  to  react  with  lithium,  and  the 
binder  reaction  could  take  on  greater  significance. 

The  simulations  elucidate  the  behavior  of  high-rate, 
lithium-ion  batteries  under  abuse  testing  based  on  the  best 
available  knowledge  of  the  reactions  that  take  place.  The 
simulations  identify  the  determining  factors  in  thermal  run¬ 
away.  Further  experimental  work  could  verify  these  factors. 
Specifically,  experimental  studies  to  determine  how  effec¬ 
tively  the  electrolyte  competes  with  the  binder  could  help 
elucidate  the  role  of  the  binder.  Biensan  et  al.  [4]  have 
reported  DSC  results  showing  the  heat  evolved  from  a 
charged  negative  electrode  is  proportional  to  the  amount 
of  PVDF  binder,  while  Maleki  et  al.  [10]  report  the  heat 
evolved  increases  with  state  of  charge  of  the  negative.  How¬ 
ever,  in  both  cases,  the  electrodes  were  dried  to  eliminate 
competing  solvent  reactions.  Repeating  these  studies  with  wet 
electrodes  (electrodes  saturated  in  electrolyte)  to  see  if  the 
total  heat  evolution  still  increased  with  the  amount  of  binder 
would  reveal  if  the  reactions  involving  PVDF  binder  can 
successfully  compete  with  those  involving  electrolyte. 

With  experimental  data,  the  parameters  in  the  simula¬ 
tions  shown  here  can  be  adjusted  to  fit  that  data.  The 


parameters  used  were  based  on  a  review  of  published 
literature,  and  oftentimes,  there  was  some  ambiguity  in 
how  to  interpret  the  parameters.  Access  to  raw  experi¬ 
mental  data  would  remove  that  ambiguity.  In  addition,  the 
simulations  carried  out  here  are  based  on  a  very  approx¬ 
imate  depiction  of  the  real  physical  system.  For  example, 
there  is  uncertainty  in  the  role  of  the  salt.  The  salt  has  been 
implicated  in  solvent  decomposition,  negative/solvent 
reactions,  and  the  positive/solvent  reactions,  but  its  role 
is  unclear.  The  mechanisms  of  all  the  reactions  are  not  well 
characterized,  and  probably  involve  several  steps,  so  the 
overall  reaction  rates  may  not  follow  simple  Arrhenius 
laws.  However,  given  the  present  understanding  of 
lithium-ion  chemistry,  the  approach  used  here  is  a  reason¬ 
able  means  to  guide  further  work. 
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